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Abstract 

This paper discusses some aspects of the global Z(N) symmetry of finite- 



temperature, SU(N), pure Yang- Mills lattice gauge theory. It contributes to 

the recent discussion of the physics of the phase of the Wilson line expectation 

ON ' value. In the high T phase, (L) takes one of N distinct values proportional to 



the Nth roots of unity in Z(N), and the Z{N) symmetry is broken. Only one 
of these is consistent with the usual interpretation (L) = e~ ' with F the 
excess free energy due to a source of color-electric flux at the position of the 



line. This relation should be generalized to (L) = ze F ' T with z £ Z{N) so 



that it is consistent with the negative or complex values. In the Hamiltonian 
description, the physical variables are the group elements on the links of the 
spatial lattice. In a Lagrangian formulation, there are also group elements 
on links in the inverse-temperature direction from which the Wilson line is 
constructed. These are unphysical, auxiliary variables introduced to enforce 
the Gauss law constraints. The following results are obtained: The generalized 
relation (L) = ze~ ' , which has appeared in earlier papers, is derived. The 
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value of z £ Z(N) is determined by the external field that is needed for taking 
the infinite-volume limit. There is a single physical, high-temperature phase, 
which is the same for all z. The global Z(N) symmetry is not physical; it acts 
as the identity on all physical states. In the Hamiltonian formulation, the 
high-temperature phase is not distinguished by physical broken symmetry. 
Rather the high-temperature phase has a percolating flux network that is not 
present in the low-temperature phase. 
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I. INTRODUCTION 

This paper discusses some aspects of the global Z(N) symmetry [|TJ of finite-temperature, 
SU(N), lattice gauge theory without matter fields. There is a low-temperature phase, which 
is presumed to confine fundamental sources, and a high-temperature phase, which does not. 
The order parameter is the finite-temperature expectation value of the Wilson line (L). The 
fundamental-representation line L(i) is the normalized trace of the product of the group 
elements on the links in the inverse-temperature direction at spatial site i. 

(L) = (L(i)) = (jjTr[UU(i,uA))) (1.1) 

The line carries a nontrivial representation of the Z(N) symmetry. In the confining phase, 
(L) = 0, and the ensemble is Z(N) symmetric. In the high T phase, (L) takes one of N 
distinct values proportional to the Nth roots of unity in Z(N). The Z(N) symmetry is 
broken. A combination of these pure phases would not have cluster decomposition for the 
Wilson line n-point functions. The broken-symmetry, pure phases appear in Monte Carlo 
simulations, where runs of finite duration in the high-temperature phase and not too close 
to the critical temperature show the Wilson line fluctuating around one of the iV available 
values. 

In a Hamiltonian approach f2|, the line is a projection operator that forces the gauge 
field to be in a fundamental rather than a singlet state at the spatial position of the line. 
Thus it is usually interpreted as the exponential of the excess free energy due to the source: 
(L) = e~ F l T . This implies that (L) should be positive. Thus, there is a problem. The 
physical interpretation should be refined so that is consistent with negative or complex 
values for (L) in the nonconfining, high-temperature phase. 

There is discussion of this point in recent papers ||. Appearing there, in various contexts, 
is a modified relation between the line and the excess free energy (L) = ze~ F ^ T with z G 
Z(N). There is also discussion of the possibility that the Z(N) symmetry is not physical 
and that the N high-temperature, pure phases associated with the N possible values for (L) 
are not physically distinct. As a consequence, well-defined interfaces may not exist. 



Here, I offer a discussion of these issues that emphasizes simplicity and gives a short, 
direct connection between basic principles and results. I use the Hamiltonian description 
of finite-temperature gauge theory on a lattice in three spatial dimensions in the Aq = 
gauge with Gauss law constraints. In this formulation, the physical interpretation is most 
accessible, and the results come fairly easily. I will refer to this as the physical theory. The 
physical variables are the group elements on the links of the spatial lattice. The constraints 
are enforced through the use of additional group elements on sites. These are unphysical, 
auxiliary variables. In a four-dimensional, Lagrangian formulation, they become the ele- 
ments on links in the inverse-temperature, 4-direction. The Wilson line is constructed from 
them. I obtain the following results: The relation (L) = ze~ F l T is verified. The value of 
z is determined by the external field that is needed for taking the infinite-volume limit in 
the high-temperature phase. There is a single physical, high-temperature phase, which is 
independent of the value of z. The global Z(N) transformation is not physical: it acts as the 
identity on all physical states. Thus, in the Hamiltonian formulation, the high-temperature 
phase is not distinguished by broken symmetry. Rather the high-temperature phase has a 
percolating flux network that is not present in the low-temperature phase, where there are 
only finite clusters of flux. 

There is considerable overlap of the present paper with Refs. |§ . Some of the statements 
in those papers are confirmed or derived by other means here. There are several additional 
points that I emphasize: It is important to use an external field and some care in the infinite- 
volume limit. The value of the factor z in (L) is determined by that external field. Even in 
a pure phase with (L) ^ 0, the local projection operators function correctly to enforce the 
Gauss law constraints. It is percolating flux rather than any physical symmetry breaking 
that distinguishes the nonconfining phase from the confining one. 

These points are discussed first in the context of a simple model for flux on links. It can 
be interpreted as a degenerate Z{2) gauge theory and is equivalent to the three-dimensional 
Ising model. 

With this Ising flux model understood, it is easy to repeat the analysis for SU(N) gauge 



theory. 

In Sec. II, the Ising flux model is discussed. Section III has the treatment of SU(N) 
gauge theory. Section IV discusses percolating flux. There is a brief conclusion in Sec. V. 

II. ISING FLUX 

This section discusses a model which is simpler than SU(N) gauge theory. However, it 
contains the essential elements of the problem at hand. It is a model for flux on links that is 
motivated by finite-temperature, SU(2), gauge theory and is equivalent to the Ising model 
in three dimensions. It can also be viewed as a Z{2) gauge theory without spatial plaquettes 
in the action. The analogue of Wilson lines on sites is spins on sites so (L) — > (s). Some 
of the critical properties of this model are discussed in Ref . 0] . Here the model serves as a 
convenient arena in which to introduce the approach that is used on SU(N) gauge theory 
in Sec. III. 

For the Ising flux model, there are results that are analogous to those for full gauge 
theories: In the high-temperature phase, (s) = ±e~ F ^ T with the sign determined by the 
external field. There is a single physical, high-temperature phase, which is independent 
of the sign of (s). The global Z(2) transformation acts trivially on physical states. The 
high-temperature phase is not distinguished by broken symmetry. Rather the high and 
low-temperature phases differ in that the high-temperature phase has a percolating flux 
network whereas the low-temperature phase has only finite clusters of flux. This last point 
was covered in some detail previously f|]. 

A. Flux model details 

In the Hamiltonian, strong-coupling description of SU(2) lattice gauge theory, the flux on 
links is labeled by the irreducible representations of SU(2). At each site, the representations 
must combine to have a singlet piece. If all of the nontrivial flux is in the fundamental 
representation of SU(2), then there must be an even number of links with flux at each site. 



The flux model results from further simplifications. First reduce the gauge group from 
SU(2) to its center Z(2). Consider it in a Hamiltonian or transfer matrix formalism. Then 
delete the contributions to the Hamiltonian from the spatial plaquettes. With that, the 
Hamiltonian is very easily diagonalized. All that remains is the energy of the links that are 
in the nontrivial rather than the trivial representation of Z(2). Those are the links with 
flux. 

In this simplified model, space is a three-dimensional cubic lattice. There are variables 9 
on links that can have the values 1 or to indicate the presence or absence of flux. These are 
the physical variables. A configuration is specified by the function 0(1). States with definite 
flux |#) are labeled by that function. A general state has a wave functional that gives the 
amplitude for the system to be found in the various basis states of definite flux. The energy 
of a link with flux is a. The weight for a configuration is 

(9\e- H / T \9) = e-iE^W = IJ e -^». (2.1) 

i 

The sum over configurations is restricted to those in which the number of links at a site that 
have flux is even. Let the collection of all such configurations be C . It is a subset of the 
unrestricted collection of all configurations C . The partition function is 

Z = Tr'[e- H ' T ] = J2( e \ e ~ H/T \ e ) = Ee^E^W. (2 .2) 

C" c 

The equivalence of this flux model to the Ising model results from using site variables to 
enforce the restriction on configurations, and then doing the 9 sums ||. The site variables 
are the Ising spins. First, consider the sum of 9(1) over the 2d = 6 Vs contained in the set 
I(i) of links with endpoint i: 

E(i) = £ 9(1). (2.3) 

To force this to be even at each site, introduce the site variables s(i), which take the values 
±1. The factor 

| E *« EW (2-4) 

s(i)=±l 



has the desired effect. A factor like this is introduced into the partition function sum for 
each site i. Now the restriction on the 9 sums can be relaxed. With the abbreviations 

E-IKIE), (2.5) 

s i s(i) 

z = e _CT//T , and ss' for the product of the pair of spins on the sites that bound a given link, 
the partition function is 

z = J2EU^ am U^) m (2-6) 

= EEiWt ( ° 1 (2-7) 

s C I 

= i; IlC 1 + *««')■ (2.8) 

8 I 

This is a version of the Ising model 

Z~^e^ ss ' (2.9) 

s 

with z = tanh/5. The links with flux are those that receive the term zss' rather than 1 in 
the expansion of Yli(l + zss'). Note that I never refer to the Ising model 1//3 as temperature. 
Temperature always refers to the T first appearing in ( |2.1|) . Small T gives small z, small /3, 
and disordered Ising spins. 

From the perspective of this paper, the physical variables are flux on links, and the 
physical states are the flux configurations. The Ising spins are unphysical, auxiliary variables 
introduced to enforce the restrictions on flux configurations. They are the analogue of Aq 
or of the gauge group elements on links in the inverse-temperature direction. 

Now consider a quantum, Hamiltonian view of the constraints. With % a site of the 
spatial lattice and g(i) G Z(2), let V(i,g(i)) be defined by its action on the basis states of 
definite flux by 

V(i,g(i))\...9(l)... >= g(i)^%..e(l)- > ■ (2-10) 

(Strictly speaking, one should distinguish whether a link contributes g or g^ 1 . However in 
Z(2), they are equal, and it does not matter.) Evidently V commutes with the Hamiltonian. 



It is the realization of a gauge transformation for this model in the definite-flux basis. A 
full gauge transformation is effected by V(g) = YliV(i,g(i)). A global Z(2) transformation 
has g(i) = g independent of i and is written V(g). 

The irreducible representation of the group Z(2)i at site i carried by \6) can be the 
trivial one corresponding to E(i) even or the nontrivial one corresponding to S(i) odd. The 
constraint on configurations (analogue of the Gauss law) says that the representation must 
be trivial at all i. Using r = or 1 to indicate the even or odd representations, respectively, 
the projection operator for the representation r(i) at i is 

P{i,r{i)) = $ y Eg r ®V{i,g). (2.11) 

9 

To select the trivial representation at all sites, use 

V = UP{i,0). (2.12) 



Comparing fl2.10| ), ( |2.11| ), and fl2.12| ) with Q2.4| ) and Q2.6| ) and identifying g with s, shows 



that the spin sums are the g sums for the projection operators. Similarly, the numerator of 
the n-point function (s(i)s(j) . . .) has projections onto the odd representation at the sites 
i,j, . . . and onto the even representation at all others, while the denominator has projections 
onto the trivial representation at all sites. This gives an interpretation of 
( S (i)s(i) ) _ Tr[e~^P^l)P(j,l)...U k ^...P(k,0)] 

as e~ F / T with F the excess free energy due to odd sites at i,j,.... 

B. Global Z(2) symmetry 

The flux model is equivalent to the familiar three-dimensional Ising model with spins 
s = ±1 on sites. There is a global Z(2) symmetry s — » — s. In the region of broken 
symmetry T > T c , there are two pure phases with (s) greater than zero in one and less than 
zero in the other. However for the flux model, the spins are unphysical so that there is not 
necessarily any physical symmetry breaking. In fact, there is not. 
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Consider this in the Hamiltonian view. The global operators V(g) act trivially on all 
states whether g = 1 or g = — 1. The first is a tautology because V(l) = 1. For the 
second, notice that in the global version of ( |2.10| ), each link appears twice and gives the 
factor g(i) e ^g(i') e ( l \ If the group elements g(i) and g(i') at the two ends of the link I are 
equal, then this factor is one. Thus, the global transformations act trivially on the states. 

That means that one can insert the operator V(— 1) into any matrix element without 
changing its value. On the other hand, 

V{-l)V{g) = V{-l)\[V{i ) g{i)) = \[V{i ) -g(i)). (2.14) 

i i 

This implies that in Z the weight for the g sums from the projectors is even in global g — > —g. 
This is the same as the Ising model global symmetry s — > — s. 

The crucial point is that V(— 1) acts as 1 on physical states so that the global Z(2) 
transformation of the spins is not a physical transformation. Thus when the spins order, 
there is no physical symmetry breaking because there is no physical symmetry to break. 

This can be viewed in another way. Think of doing the 9 sums first with fixed spins (or 
g's) and then doing the spin sums. The weights for the flux configurations depend upon the 
values of the spins. However, if all the spins are reversed, there is no effect on the weights of 
the flux configurations in ( p.7|) . When the global Z(2) is broken, there are two pure phases 
for the spins, but the configurations of flux are independent of which is selected. Thus the 
physical variables are not affected by the selection of a single pure, spin phase as is needed 
for cluster decomposition in spin variables. This is the analogue of the statement that the 
gauge theory pure phases with different values of (L) are physically indistinguishable. In 
calculating physical matrix elements, one could just as well sum over the two pure phases. 
This would give the same physical matrix elements but (s) = 0. 

So if it is not physical symmetry breaking that distinguishes the confining and the non- 
confining phases, then what is it? It is the presence of a percolating network of flux. The 
physical percolating network allows the nonphysical spins to order. This will be discussed 
further in Sec. IV. 



C. Physical interpretations 

Since the physical states of flux are the same for either sign of (s), how should the sign 
and magnitude of (s) be interpreted? As shown, the sum on the Ising spin at a site can 
be interpreted as projection onto the subspace of states with even flux at that site. In the 
same way, the expectation value (s(i)) or the correlation function (s(i)s(j)) projects onto 
states with odd flux at the indicated sites and even flux elsewhere. As shown in ( |2.13| ), these 
objects are sums of diagonal matrix elements of projection operators. They are interpreted 
as e~ F l T where F is the excess free energy of the ensemble with odd flux at some sites. Thus, 
these expectation values should be positive. This is in conflict with the fact that (s) can be 
less than zero for T > T c . The apparent contradiction casts some doubt on the action of 
the supposed projection operators when the spins order. Evidently the problem is related 
to the infinite-volume limit. In finite volume, there is no ordering of the spins. 

One approach is to finesse the problem by ignoring the single spin expectation value 
and looking at the two-point function, which is even in the reflection operation and does not 
suffer the problem of interpretation. The two-point function (s(O)s(i)) has the same value in 
the two broken-symmetry, pure phases. The interpretation (s(O)s(i)) = e~ F2 ^ l '' T with ^(z) 
the excess free energy of two odd sites separated by distance i is consistent. Take i —* oo 
to obtain F 2 (oo) = IF and e~ F ' T = (s(O)s(oo)) 1 / 2 = ((s)(s)) 1/2 = \(s)\. This confirms 
the relation between F and \{s)\. However, it leaves the sign of (s) uninterpreted and the 
contradiction unresolved. 

A technique for handling this kind of problem is to introduce an external field h and take 
limits in the order h — > after volume V — > oo. The trivial action of V(— 1) on physical 
states gives the reflection symmetry for the spins. It follows that the free energy is an even 
function of h. When the volume is finite, there is an expansion in powers of h, (s) — > as 
h — > 0, and there is no problem. However for T > T c and V — > oo, the free energy has a 
nonanalyticity of the form \h\. This gives (s) ~ e(h). It shows that the sign of (s) simply 
reflects whether h — > + or h — >• 0~. 
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It was shown above that the configurations of flux are unaffected by a global reversal 
s — ► —s. It follows that the weights for flux configurations are unaffected by h — ► — /i. Thus 
whatever (s) has to say about flux, it says the same thing for h — ► + or /i — ► 0~. 

So what does (s) say about flux? For h = in finite volume, (s) = 0, and the interpre- 
tation is very simple: Projection onto a state with one odd-flux site is zero because it is a 
geometrical impossibility to have exactly one odd-flux site. The h — 0, infinite-volume limit 
does not exist. If it did, it would surely give (s) = 0, which is not correct if spin cluster 
decomposition is desired. With h ^ 0, the infinite-volume limit is defined. The external 
field gives an extra factor of e h ^i s ^'. To interpret it, write out the character expansion 

e *£i«W = JJ e h»W (2.15) 

= Y[(coshh + s{i) sinh/i). (2.16) 

i 

This shows that the external field inserts odd-flux sites with a density proportional to tanh h, 
which vanishes as h — ► 0. 

In finite volume, there is a small-h expansion that begins 

(s(0)) h = tanh/i£<s(0)s(z)) + ... . (2.17) 

i 

In this expansion, the n-point functions on the right hand side are evaluated at h — 0, and 
only those with even n appear. For odd n, the correlation functions in finite volume with 
h = are zero. 

The leading, large- volume terms sum to a function of the form t(Vh). The function t(x) 
must have a small x expansion with odd powers of x and a large x limit with £(±oo) = ±(s). 

This shows that (s) is determined by the infinite-volume limit of n-point functions with 
n even. These are unaffected by global sign changes so the sign of (s) is entirely determined 
by the sign of h. Since the flux configurations are the same for either sign of h, the sign of 
(s) is unphysical. 

When the spins order, one pure phase with a definite sign for (s) is selected. One does 
not ordinarily sum over the two pure phases. Does this mean that some global part of the 
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constraints or projections is not being handled correctly? No, because the action of V(—l) 
on physical states is trivial. 

For further peace of mind, one can verify that the local projectors are functioning prop- 
erly. This means that there should be contributions only from configurations with even flux 
at each site in a single, ordered, pure phase. A measure of this comes from the operator 
J-(i) = V(i, —1) = (— 1) E W. If the expectation value of this is 1, then only even-flux sites 
are contributing. Inserting this operator is the same as reversing the spin at site i. Thus, 

(Hi)) = git^k^p^ . (2.18) 

But a change of summation variable s(i) — *■ —s(i) on the top reduces this ratio to 1 whatever 
the phase. For more detail, see Sec. III.C. 

D. Lagrangian view 

For a full analogy with ordinary gauge theory, this Hamiltonian formulation in three 
spatial dimensions can be related to a 4-dimensional Lagrangian formulation in which the 
lattice has finite extent and periodic boundary conditions in the inverse-temperature, 4- 
direction. Consider a gauge theory with elements of Z(2) on the links of this 4-dimensional 
lattice. Choose an action with the usual contributions from the plaquettes with links in 
the inverse-temperature direction but with no terms from spatial plaquettes. Now choose 
a gauge where all 4-direction link variables are 1 except those on a particular slice at fixed 
i&. These are the Ising spins. To obtain the effective theory for these spins, sum over the 
variables on spatial links. An Ising interaction among the spins remains. Thus, the Ising 
model is a Z(2) gauge theory without terms in the action from spatial plaquettes. 

In this Lagrangian view, there is a global Z{2) symmetry operation that multiplies all 
the variables on 4-direction links at a given i 4 by — 1. Viewed as a classical theory on the 
4-dimensional lattice, the finite-temperature gauge theory has a global Z(2) symmetry. 

This is the most common description of finite-temperature gauge theories. Viewed as 
classical theories, there is a Z(2) symmetry that can be broken. But, as discussed above, 
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when they are viewed as quantum systems in three spatial dimensions at finite temperature, 
there is no physical Z{2) symmetry. 

E. Ising flux conclusion 

For the Ising flux model, whatever the phase, only configurations with even flux at each 
site contribute to the sums for the partition function. The weights of the flux configurations 
are the same in the two pure phases of ordered spins. There is no physical difference between 
these phases. The expectation value of a spin is (s) = e(h)e~ F ^ T . The quantity F is the 
excess free energy for the ensemble with a source at one point and a background of other 
sources in the limit of vanishing density for the other sources. The sign, which is determined 
by h, is an artifact of the limiting process and has no physical significance. 

III. SU(N) GAUGE THEORY 

Now that the concepts have been developed in an elementary setting, it is relatively 
easy to generalize to SU(N) gauge theory. First there is a short paragraph on the usual 
Lagrangian description of the finite-temperature, Z(N) symmetry and its breaking. Then 
there is a discussion using the Hamiltonian machinery. That is followed by a brief mention 
of the analogue of the Ising spins view. 

The main difference from the Ising flux model is that states of definite flux are no longer 
eigenstates of the Hamiltonian or of the transfer matrix. However, that is not an essential 
part of the argument. The important thing is that they are a basis, and that they still are. In 
fact, it is only a convenience to work in that basis. The arguments can be made just as well in 
the other commonly used basis where the group elements on links have definite values. The 
global transformation by g — — 1 acts trivially on physical states, but it produces the global 
Z(N) symmetry of finite-temperature gauge theory. It is a symmetry of the unphysical, 
auxiliary group elements on links in the 4-direction. So it is not physical symmetry breaking 
that distinguishes the confining from the nonconfining phases. This is evident in either basis. 
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However, the definite-flux basis is clearly more convenient in describing the percolating flux 
that does distinguish the confining from the nonconfining phases. 

In the Lagrangian view, finite-temperature gauge theory is a classical system of fields 
on a 4-dimensional lattice with one finite, periodic direction. The group elements on the 
4-direction links have no special interpretation. There is a global symmetry that results from 
multiplying all the group elements on links in the 4-direction at a given i& by an element 
of the center of the group. The breaking of this symmetry, as detected by the Wilson 
line, reveals the deconfining phase transition. To give this a quantum, finite-temperature, 
physical interpretation, one requires a Hamiltonian description. 

A. Hamiltonian Machinery 

The Hamiltonian description of the finite-temperature gauge theory on a lattice in 
three spatial dimensions in the A$ = gauge with Gauss law constraints is what I refer to 
as the physical theory. The physical variables are the group elements on the links of the 
spatial lattice. The basis states for the quantum system can be taken as either those with 
definite flux, which is definite representation functions on links, or those with definite group 
elements. 

This state space carries a representation of the gauge group G = l\i SU(N)i. As before, 
let V(i,g(i)) be the operator for a gauge transformation by g(i) G SU(N)i at site i. To 
enforce the Gauss law constraint at site i, one projects onto the subspace that carries the 
trivial representation (now indicated by 1) of this group action. The operator that does this 
is 

P(i,l) = Jdg(i)V(i,g(i)). (3.1) 

If there are no external sources, then V — Yli P{h 1) projects onto allowed states. The 
partition function is 

Z = Tr[e- H/T V]. (3.2) 
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When Z is converted to a path integral, the g(i) integrals from (|3.1|) become the integrals 
over the group elements on the 4-direction links. 

To account for an external source at site i in representation r(i), the projector P(i, 1) in 
the product V is replaced by one with a factor of the character for irreducible representation 



r[i) 



P(i, r(i)) = Jdg(t) x*r {i) (9(i)) V(i, g(i)). (3.3) 

In the path integral, this becomes a Wilson line in representation r(i) at site i. In this paper, 
only the trivial representation r = 1, the defining, fundamental representation r = N, and its 
complex conjugate r = N* are needed. All sites are projected onto the trivial representation 
except at most a few onto these fundamental representations. The line expectation values 
are expressed in terms sums of diagonal matrix elements of projection operators. 

mm ■) = ^»^^ p., 

In this view, the physical variables that describe a state of the quantum system are the 
group elements on spatial links. The group elements on 4-direction links are unphysical, 
auxiliary variables used in a technical way to enforce the Gauss law constraints. 

B. Global symmetry 

Now let z stand for one of the N elements of the center Z(N) of SU(N). Consider the 
global transformation V(z) = YliV(i,z). In either basis, this is easily seen to act trivially 
on the physical variables. For example, with definite group elements u(l) on spatial links I 
bounded by sites i and i', 

u(l) -► g- x {i)u{l)g{i') -► z^utyz = u(l). (3.5) 

Similarly, when the flux is definitely in representation r(i) with indices m and n, the am- 
plitude for the group element on link I to have a value u(l) is given by the representation 
matrix D^(u(l),r(i)). The relations 
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D(z- X u(l)z, r(i)) = D-\z, r(i))D(u(l),r(i))D(z, r(i)) = D(u(l),r(i)) (3.6) 

show again that the action of V(z) is trivial. In the absence of a nontrivial, physical, global 
transformation, there can be no physical symmetry breaking at the phase transition. The 
confining and nonconfining phases are not distinguished by different physical symmetry 
realizations. 

Since the operator V(z) acts as 1, it can be inserted into any matrix element without 
changing its value. However, because 

V(z)V(g) = V(z) I] V(i, g(i)) = U V(i, zg(i)), (3.7) 

i i 

there is the effect of changing each auxiliary variable from g(i) to zg{%). With fixed g's, the 
weight for a flux configuration is 

(e\e- H / T UV(i,g(i))\0)- (3-8) 

i 

Under the global transformation g — > zg, it is unchanged: 
(9\e- H / T l[V(i,g(i))\e) - (e\e- H / T l[V(i,zg(i))\0) = (0\e- H / T V(z) l[V(i,g(i))\9) (3.9) 

i i i 

= (0\e- H/T UV(i,g®)\0)- (3-10) 

i 

In the path integral, this becomes the global Z(N) symmetry of the 4-direction link variables 
in the classical, Lagrangian view. Thus, the ordering of the Wilson lines does not break any 
physical symmetry. 

Furthermore, a similar argument using V(z) gives 

o = JdgxUg)(<P\v(i,gM) (3.n) 

for any physical states and ip. Thus, any calculation that gives (L) ^ cannot be 
interpreted as 

Tr[e- H / T P(i,N)HP(j,l)}. (3.12) 

To find the correct interpretation of (L) ^ 0, the external field is introduced again. 

16 



C. Interpretations 

At this point, the relation of the global symmetry of the Wilson lines to the physical, 
Hamiltonian formulation has been given. As indicated in ( |3.10| ), when the auxiliary variables 
g(i) that enforce the Gauss law constraint and become the Wilson lines are multiplied by 
an element of the center of the group, the weights for the different physical configurations 
of flux are unchanged. Thus, the physical configurations are insensitive to the phase of (L) 
in the broken-symmetry region. Each of the N ordered, pure phases gives the same result 
for physical quantities. In calculating physical quantities, one could just as well sum over 
these pure phases and have (L) = 0. 

To treat this more carefully and find an interpretation for (L) ^ 0, an external field 
h is introduced again. Let L(i) be the line at site % computed in the defining, funda- 
mental representation N. The logarithm of the Boltzmann weight gets the additional term 
h*J2iL(i) + hJ2iL*(i). Now g — > zg is equivalent to h — > z*h, so the free energy must 
be invariant under a rotation of h by an element of Z(N). Conversely, modifying h by 
such a factor in the Boltzmann weight has no effect on the weights of the different physical 
configurations. 

In a finite volume with h — > 0, there is no singularity, and the Z(N) invariance guarantees 
(L) = 0. In that case, there is no contradiction to resolve. The problem is in infinite volume 
where (L) ^ is possible. An interpretation of (L) as the expectation value of a projection 
operator or as e~ F / T requires (L) > 0. Also, Z(N) invariance requires (L) = 0. Thus, a 
negative or complex value for (L) requires some explaining. 

Of course, the problem can be avoided again by ignoring it and extracting the desired 
information from the two-point function at infinite separation (L(0)L*(oo)). But this does 
not do much to refine our understanding or resolve the contradiction. 

With an extra factor of e h ^» L W +h ^ L W inserted in the Boltzmann weight, there is 
a finite density of sites projected onto nontrivial representations. Again, the character 
expansion makes this precise. At each site, the factor is 
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Ji* L{i)+hL* (i) 



"£a r (h,h*) Xr (g(i)) (3.13) 

r 

OlU + E-Xr) ( 3 - 14 ) 

ai(l + — Xn + — X^) + 0(h 3 ) (3.15) 



with 



h*/N + 0(h 2 ) and — = V^ + 0(/i 2 ). (3.16) 

a x ai 

The density of sites in the fundamental representation is 0(h*). 

Thus, the expectation value of the line is computed in the presence of a finite density of 
other sites in nontrivial representations. As h — > 0, this density goes to zero. The small h 
expansion starts as 

(L(0)) h = hJ2(L(0)L*(i))o + ... • (3.17) 

i 

The large- volume limit of this function depends on the phase of h so that as h — > 0, the 
result is proportional to the nearest element of Z(N). For example, 

{m)h = ^W lnZ (3 - 18) 

with 

Z= Y. e v ^ {h *" z+hz '\ (3.19) 

zez(N) 

The large-volume limit depends on the phase of h. However, each matrix element that 

goes into that limit of Q3.171 ) is computed in finite volume with no external field and is a 

global Z(N) singlet. The density of the nontrivial sites goes to zero, and those with i ^ 

recede to infinity as h — > 0. The contribution from each flux configuration is unchanged 

by the replacement h — > zh. Thus, the flux configurations and all physical properties are 

independent of the phase of (L) . The conclusion is that 

(L(i)) = ze~ F/T . (3.20) 



F is the excess free energy associated with a source at i in a background of other sources as 
their density goes to zero. The phase z G Z(N) is determined by the direction from which 
the external field approaches zero. 

To confirm that there is no change in the realization of the gauge properties of the theory 
in the nonconfining phase, one needs to check that even though the global Z(N) is broken, 
there is no effect on the local Gauss law constraints. Since it has already been shown that 
global V(z) has no physical effect, this is a pretty safe bet. 

The constraint is supposed to project onto the trivial representation of SU(N)i. If that 
is true, then 

(V(i,g')) = l . (3.21) 

The effect of this extra operator is to multiply the auxiliary variable g(i) by g'. However, 
that effect can be immediately removed by a change of variable in the g(i) integration. Thus, 
( P-21 ) is correct, and the Gauss law is implemented. This works whether or not the g's have 



ordered. 

If you are nervous about that in the ordered phase, consider it more carefully in the 
presence of the external field. In that case, there is additional g(i) dependence in the weight 
from e h ^* w+ h 2^n L W. Since L(i) is the trace of g(i) in the defining representation, when 
the integration variable is changed, there is a new factor 

e h*{z-l)L{i)+h(z*-l)L*{i) (3^2) 

associated with site i only. This is a local change, so the global phase is still determined 
by the field h acting over the whole volume. The factor (|3.22|) becomes 1 as h — > 0, the 
expectation value of V is 1, and the local Gauss law is correctly implemented even in an 
ordered, pure phase. 

To complete the parallel with the Ising flux model, the analogue of the Ising spins view 
must be mentioned. For gauge theory, this view is known as the effective theory for the 
Wilson lines || . Holding the values of the lines fixed, all the other variables are integrated 
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out. This leaves a theory in three space dimensions with fields L{i) on sites. There is a 
global Z(N) symmetry. 

D. Monte Carlo 

Monte Carlo calculations are ordinarily done in finite volume with no external field. How 
should they be interpreted? If one were to run long enough, the result would be (L) = 0. 
However, in practice, in the nonconfining phase and not too close to the phase transition, the 
flips of L from one pure phase to another occur infrequently or not at all so that (L) ^ for 
a run. This has no effect on the measurement of physical variables. A series of configurations 
in which the line does not flip corresponds to the h — ► after V — ► oo procedure. They 
must be interpreted according to the discussion in the previous subsection. They cannot be 
given the simple interpretation 

(L(i)> ~ Tr[e- H l T P{i, N) \{ P(j, 1)]. (3.23) 

IV. PERCOLATING FLUX 

I have asserted that the difference between the confining and the nonconfining phases is 
not a difference in the realization of a physical symmetry but rather the presence or absence 
of a network of percolating flux. Let me elaborate. 

A. Ising flux 

It is helpful to look at Ising flux first. Consider two widely separated points % and j 
and the correlation function (s(i)s(j)). In the low-temperature phase, flux is sparse, and 
the character expansion is convergent. If the sites i and j are not connected by any flux, 
then the spin sums in the separate clusters that contain % and j are independent, and the 
contribution to (s(i)s(j)) is zero. When the probability of arbitrarily large flux clusters 
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goes to zero as their size increases, the correlation function must vanish as the separation 
becomes arbitrarily large. Thus, no percolating cluster implies no ordering. 
On the other hand, if there is ordering so that 

lim (s(i)s(j)) = (s) 2 > 0, (4.1) 

\t-3\-*oo 

then there must be a finite probability to find i and j in the same cluster no matter how 
large \i — j\, i.e. there is a percolating cluster. Spin ordering implies a percolating cluster. 

Conversely, if there is a percolating cluster of finite density, then there is a finite prob- 
ability to find both % and j in it no matter how large \i — j\. It is necessary to consider 
the \i — j\ dependence of this probability. I do not know how to address this rigorously, 
but consider the following discussion: In any finite volume, the largest cluster hits a finite 
fraction of the sites. If that fraction remains finite as the volume goes to infinity, then 
there is a percolating cluster. It seems clear that energy and entropy will conspire to give 
the percolating cluster an approximately constant density. Large departures from constant 
density have low entropy. With energy and entropy conspiring to give constant density for 
the percolating cluster, the disturbance caused by an odd-flux site will extend over a fixed, 
finite distance. If i and j are separated by more than that, then further separation will have 
no effect. Assume that the weight for a largest cluster is a function of its density. Since 
there are largest clusters with the same density independent of \i — j\ when it is large, the 
probability to find both % and j in the largest cluster is also independent of \i — j\. 

After the spin sums at sites other than % and j have been done, one is left with a factor 
proportional to [s(i)s(j)] as the contribution from all the links in the largest cluster. The 
final sums are 

(^)#~EE^)[#(i)]^)>0' ( 4 - 2 ) 

s(i) s(j) 

Finally, take the infinite- volume limit. This shows that percolating flux implies spin ordering. 
Note that the choice of a pure phase for the spins causes no difficulties because the two-point 
function and the flux distributions are independent of that choice. 
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In conclusion, the physical difference between the nonconfining and the confining phases 
is the presence of percolating flux. 

B. Gauge flux 

After the discussion above, it is easy to apply the ideas to gauge flux. The flux on links is 
more complicated since it is labeled by the irreducible representations of SU(N) rather than 
Z(2). However, the argument above is geometrical in its essentials and does not depend upon 
the detailed nature of the flux. It is necessary that the auxiliary variable group elements g(i) 
and g(j) be uncorrelated if % and j are in different clusters. This is certainly true. Define 
"there is a percolating cluster" to mean that percolating clusters make a finite contribution 
to (L(i)L(j)) when the separation is finite. Now, no matter how large the separation, there 
is a finite probability to have both sites hit by the percolating cluster in a flux configuration. 
If the disturbance in the cluster caused by one nontrivial site has a finite range, then the 
value of the contribution from the percolating clusters is independent of \i — j\ for large 
\i — j\. Thus, the reasoning from the Ising flux case can be used again. The conclusion is 
that the physical difference between the nonconfining and confining phases is the presence 
of percolating flux. 

V. SUMMARY 

The usual relation (L) = e~ F ^ T cannot be correct for nonzero, nonpositive values of 
(L). A careful approach to the Z(N) symmetry and to the infinite- volume limit resolves the 
problem. In terms of the physical variables of quantum, finite-temperature, Hamiltonian 
gauge fields, there is no Z(N) symmetry to be broken. The physical system is the same in 
each of the N pure phases that can be chosen when (L) ^ 0. The confining and nonconfining 
phases are distinguished by the presence of percolating flux in the nonconfining phase. To 
treat the infinite-volume, (L) ^ pure phases, an external field is introduced in the usual 
manner. This leads to the conclusion that (L(i)) = ze~ F / T . The element z G Z(N) is 

22 



determined by the external field. F is the zero-density limit of the excess free energy of a 
fundamental source of flux at site % in a background of other sources. 
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